Parametrization of spectral surfaces of a class of 
periodic 5-diagonal matrices 

Ionela Moale and Peter Yuditskii* 
December 15, 2011 

Abstract 

The main result of this work is a parametric description of the 
spectral surfaces of a class of periodic 5-diagonal matrices, related to 
the strong moment problem. This class is a self-adjoint twin of the 
class of CMV matrices. Jointly they form the simplest possible classes 
of 5-diagonal matrices. 
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1 Introduction 

The main result of this work is a parametric description of the spectral 
surfaces of a class of periodic 5-diagonal matrices. This class is a self-adjoint 
twin of the quite popular now class of CMV matrices [5, 6]. Jointly they 
form the simplest possible classes of 5-diagonal matrices, we explain this 
claim in details in Section 2. We take a risk to introduce an abbreviation 
and call them SMP (Strong Moment Problem) matrices, see e.g. [3] and the 
remark below. 

In a generic case SMP matrix is defined as a ratio of two 3-diagonal 
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matrices A = A2A l 1 , where 
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A\ and A 2 are invertible, equivalently &m-\ 7^ and q<m 7^ for all n. Due 
to the symmetry A = A*, the generating coefficient sequences are subject 
to the restrictions 
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Note that Aj" 1 and 1 are also 3-diagonal, thus A and A -1 are 5- 
diagonal matrices represented as products of 3-diagonal matrices. Note also 
that degenerations are possible, see Remark 4.3. 

Let {e„} be the standard basis in the two-sided £ 2 . According to the 
above definition 

Ae2n = Aiein = P2n^2n-l+<i2ne-2n + P2n+l e 2n+l 
A 1 e2 n -l = A\e2n-\ = ^2n-l&2n-2+^2n-l ( i2n-l + ^2n,e2n 

Thus, constructing A we follow the procedure, which is similar to the CMV 
matrices case: having e_i and eo as the generators of the cyclic subspace we 
form the whole space applying A on the even step and A^ 1 on the odd step, 
however, as it was mentioned, the operator is not unitary, but self-adjoint, 
that is, the spectrum is not on the unit circle but on the real axis. From this 
point of view the situation is similar to the orthogonalization procedure in 
the strong moment problem construction [3] , however because of periodicity 
we are interested in two-sided matrices; it is very essential: we do not assume 
that the operator A is positive! 

A can be represented as a two dimensional perturbation of a block or- 
thogonal matrix 
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where 



po = ||-P + Ae_i || , e = —P+Ae-i, 

Po 



A± = P±AP± are restrictions of A to the positive and negative half- axis 
according to the orthogonal decomposition I 2 = l 2 _ © t\. For this reason 
certain general facts from its spectral theory can be reduced to the spectral 
theory of Jacobi matrices. However it is quite different as soon as we pose 
the problem: 

Problem 1.1. Describe the spectral sets of periodic SMP matrices. 

Recall that for periodic Jacobi matrices the spectrum is a system of 
intervals, which possesses the following parametric description, see [1] and 
references therein. For a system of nonnegative parameters {/ifc}£zj, let 
D = D(hi, . . . h n -i) be the region obtained from the half-strip 



Let 9 be the conformal map of the upper half-plane H to D normalized by the 
conditions 9(ao) = 0, 9(bo) = —irn, 9(oo) = oo. Denote by E(h±, . . . , h n -\) 
the full preimage of the interval [— 7m, 0] C dD, i.e.: 



A system of intervals E = [bo, ao] \ L)j>i(aj, bj) is the spectrum of a periodic 
Jacobi matrix if and only if E = E(h\, h n -\) for a certain system of 
parameters {/&fc}£zj. 

The spectral sets for periodic CMV matrices are given by conformal 
mappings onto similar periodic comb-like domains. To formulate our main 
result we define comb regions of a new kind. 

For integers k and m G [0, k] and parameters ho > and ojq, < ujq < 
ixm consider the region n™(/io,^o) bounded by the hyperbolic curves 



{w : —ttti < Rero < 0, Imw > 0} 



by removing vertical intervals 



{w : Rew 



irk, 0<lmw<hk}, k = 1, . . . , n — 1. 



E(h 1 ,...,h n - 1 ) :=9- 1 ([-nn,0}). 



(Imu) 2 < (Reu) 2 + h 2 , 



(1.1) 



and the orthogonal systems of hyperbolas 




for Reu < 0. 



for Ren > 0, 



(1.2) 
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Figure 1: IT region for luq / tt£ 




Figure 2: II region for loq = tt£q. 



4 



If ojo / ir£, < £ < m, by II we denote the region which is obtained from 
Il™(/io, ujq) by removing pieces of hyperbolic curves 

Imn = — — — , for Imu > 0, 1 < £ < m — 1 (1.3) 
2Reu ~ ~ v ' 

and 

Imn = ?r ^~ ^ for Imu < 0, m+ 1 < I < 2k- 1 (1.4) 

2Reu 

of length /i^, /i; > 0, see Fig. 1. If ujq = tt£q then the hyperbolic curves 
related to £ = £o in (1.3) and £ = k + £o in (1.4) degenerate. In this case the 
corresponding cuts are pieces of the imaginary axis, as soon as 

h to < h and h k+to < h . (1.5) 

Otherwise one of them still satisfies (1.5), and another one has T-shape, see 
Fig. 2, consisting of the piece of the imaginary axis 

< Imn < ho or — ho < Imu < 0, (1-6) 

respectively, and of the real interval 

-h { !° ] < Ren < hf\ h { l o) > 0. (1.7) 

Theorem 1.2. For a region II described by the conditions (1.1) —(1.7) , let 
6 : H — > II be a conformal map such that and oo correspond to the infinite 
points in II. A system of intervals E, 0, oo E, is the spectral set of a 
periodic matrix of SMP class if and only if it corresponds to the preimage 
of the part of the boundary given by the condition (1.1) for a certain II. 

The structure of the work is as follows: the simplest possible spectral 
surfaces of periodic 5-diagonal matrices are discussed in Section 2. In Section 
3 we prove our main theorem. The functional model for periodic SMP 
matrices is given in Section 4. 



The second author is thankful to the organizers of the program Hilbert 
spaces of entire functions and spectral theory of self-adjoint differential op- 
erators, at CRM, Barcelona, 2011, and to Alex Eremenko for stimulating 
discussions. In a sense this paper is an addition to their joint work [1]. 
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2 Spectral surfaces with the maximal number of 
boundary ovals 



Let J be a 5-diagonal self- adjoint matrix of period d 
J = rS 2 + pS + q + S~ l p + S- 2 r, 



(2.1) 



where S is the shift operator and p, r, q are diagonal matrices of period d, 
such that r m > 0. We recall certain fundamental facts from the spectral 
theory of multi-diagonal periodic matrices [2], adopting to the 5-diagonal 
case. 
For 
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(2.2) 



let 



F(z,w) = 



det{ j(w) — z ■ 1} 



n 



d-l 



w 2 + 1/w 2 + A(z)w + A*(z)l/w + B(z) 



where A and B are polynomials, in particular, for even d = 2k 

„2k 

B(z) 



A*{z) := A(z) 



( -1 



(2.3) 



+ 



nk—l T~fk— 1 

j=0 r 2j lli=0 r 2i+l, 



z k + 



Then the spectral curve corresponding to J is of the form 

U = {P= (z, w) : F(z, w) = 0}. (2.4) 
1Z is endowed with an antiholomorphic involution tP := (z, 1/w) for which 

TZ\dTZ + = TZ + = {P = (z,w) G 11 : \w\ < 1}, 

see Fig. 3. Note that the spectrum of J (as the operator acting in £ 2 ) 
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Figure 3: Topology of the spectral curve 



corresponds to the fixed line of the involution r, tP = P, 



z G a (J) <S> 3w : P = (z, w) G &R+. 



In other words it is described by the condition = 1. 

Recall that the spectral surfaces related to periodic Jacobi matrices are 
of the form 



where A is a real polynomial. In the similar decomposition 1Z \ cflZ = 
it possesses the following property: the number of boundary ovals, 
i.e., the number of intervals 



is maximal for the given genus of the surface. 

We say that the spectral curve 1Z, related to a 5-diagonal matrix, is of 
the simplest structure if it has maximal possible number of components of 
the boundary for the given genus. For example, in Fig. 3 the number 
of boundary components is 3, but its genus is 4 and the maximal possible 



1Z = < (w, z) : w + 



w 
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dn = {z£R: \A(z)\ < 2} 
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Figure 4: 1Z + as the two sheeted covering of z-plane 

number of components is 5. That is, the curve of this structure does not 
belong to the class. 

In other words, let us represent 1Z+ as two a sheeted covering of the z- 
plane, see Fig. 4. It is a hyperelliptic curve with a system of cuts &R + . We 
say that the spectral curve is of the simplest structure if this hyperelliptic 
curve has genus 0, i.e.: 

TZ + ~ C \ E. 

The corresponding equivalence can be written explicitly 

c 2 

z = A + — , A G C, ci = A - 2c, c 2 = A + 2c, (2.5) 

A - Ao 

where ci, c 2 denote the only two possible critical values of z, in the case 
that both numbers are finite and if, say, C2 = oo then 

z = A 2 - 2c, ci = c 2 - 2c. 

The set E, which corresponds to dlZ + , is a system of cuts in the complex 
plane C with the property 

z(X) G R for A G E. (2.6) 

It is essential to note that E is far from being an arbitrary system of 
cuts for which (2.6) holds. Recall that up to now the second function w was 
not involved into considerations. Meanwhile w = w(X) is a function in C\E 
with the following properties [2]: 



8 



(i) w is single-valued and holomorphic, 

(ii) \w\ < 1 in C \ E and = 1 on E, 

(iii) zeros of w are {Ao, 00} = z~ 1 (oo) <£_ E (of equal multiplicity). 

For defmiteness, here and below, we consider the case (2.5) (with two finite 
critical values). The properties (i)-(iii) imply that 

\ log r4u\=Gx (\) + Goo(A), (2.7) 

where G\ (A) is the Green function in the domain C\E with a logarithmic 
pole at Ao and k is the multiplicity of w in Ao and 00 respectively. 

Let us recall the concept of the complex Green function b\ (\). It is an 
analytic multivalued function in C \ E such that 

log sk = GA ° (A) - <2 ' 8) 

Note that (2.8) determines b\ up to a unimodular constant. In what follows 
we assume the normalizations 6a (oo) > and 600 (Ao) > 0. 

Let 7Ti (C \ E) be the fundamental group of this domain. Then b\ gen- 
erates the character n\ € vri(C \ E)* on this group by 

K°7 = i"Ao(tK. 76 7Ti(C\.E7), 

which indicates the multivalued structure of the complex Green function. 
Moreover, let us split E into connected components, E = Uj =0 Ej, and let 
7-, 's be simple contours around E/s. Note that they form generators of the 
group 7Ti(C \ E) subject to the condition 

7o • • • 7m = trivial. 
Then, for a suitable choice of the direction of jj , 

^Ao(7 J ) = e 27r ^o(^), (2.9) 

where ui\ {Ej) is the harmonic measure of Ej at Ao- 

The factor b\ boo removes the singularities of z in C \ E. Let 

(boob\ z)(oo) (boobxgz)^) 



1(6006^^(00)!' Ao |(6 00 6 Ao z)(A )|' 

Define £ £ [0,1) by the condition = e 2m ^t\ . Let [i = /J.\ fJ,oo- Thus 
A» (7j) = e 2 ™J, where Uj = uj Xo (Ej) + oo^Ej). 
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Theorem 2.1. Let z be given by (2.5). Let E = Uf =0 Ej C z _1 (K) be a 
system of cuts (2.6). Then K = TZ+ U dTZ + U where 1l + ~C\E, is 
a spectral curve of a 5- diagonal periodic matrix if and only if the numbers 
£ and ujj (for all j) are rational. Moreover w = t^ : (b\ b oo ) k , where k is a 
common denominator of these rational numbers. 

A proof is based on a fact of the general theory [2] , that a given periodic 
J with the spectral curve 1Z possesses functional representation as the mul- 
tiplication operator by z. We give such a representation following basically 
[4,7]. 

For a fixed character a the multivalued analytic functions F, F o j = 
a{-f)F, such that |F(A)| 2 has a harmonic majorant in C\E, form the Hardy 
space H 2 (a) C L 2 ^ with the norm given by the integral of the boundary 
values: 

||F|| 2 = [ |F(A)| 2 Woo (dA). 

JE 

Note that the point-evaluation functional is bounded in this space and there- 
fore in H 2 (a) there is the reproducing kernel fc": 

F(\) = {F,k%), \£C\E, 

for all F G H 2 (a). 

Lemma 2.2. For a character a G 7Ti(C \ E)* let 

denote the normalized reproducing kernel at A. Then, for an arbitrary system 
of unimodular constants t rn , the family 

— n —n — n — 1 —n 

e2n - t 2n b Xo b 00 K Xo , e 2n +i - hn+ib Xo b^Koo (2.10) 

forms an orthonormal basis in H 2 (a), n > 0. Moreover, extended to negative 
indexes it forms an orthonormal basis in L 2 . 

Proof. The system is orthonormal. Every function from H 2 (a) orthogonal 
to it has a zero of infinite multiplicity in Ao (and oo) and therefore vanishes 
identically. To prove the second claim one has to use the description of the 
orthogonal complement L 2 ^ © H 2 (a) by means of the Hardy space [7] and 
once again apply the same argument related to the corresponding .ff 2 -space 
and orthonormal basis of reproducing kernels in it. 

□ 
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Lemma 2.3. The multiplication by z with respect to the basis (2.10) is a 
5-diagonal self-adjoint matrix, 

Ze m +2 = r m +2em+Prn+2e m +l+<lm+2e m +2+Prn+3e m+ 3 + r rn+ 4e m+ 4 : . (2.11) 

Moreover r m > if and only if 

*2n-l = *^o n *(-l)> t 2 n = t^t( y (2.12) 

Proof. Since the factor b\ boo removes the singularities of z in C \ E the 
function zb^b^F belongs to H 2 (aii) for every F £ H 2 (a) and an arbitrary 
a G 7Ti(C \ -E)*. Thus the decomposition of ze m+ 2 starts with e m , 

Z&m+2 — fm+2^m ~\~ ■ ■ ■ 

Since z(X) is real on E the multiplication operator is self-adjoint; its ma- 
trix possesses the symmetry property and therefore it is 5-diagonal (2.11). 
Finally we put A = Ao in (2.11) for even m 

(^K) 6 oo)(A )^" ( f ( + ) (A )t2n+2 = f 2n +2Kxo ( A o)*2n 

and A = oo for odd m 

(z6 A0 6 oo )(oo)^^ ( " +2) (oo)t 2n+1 = f 2 ^.i^S t «" t ""(oo)t 2n _i. 

Since K<*(\) > we get (2.12). □ 

Proof of Theorem 2.1. Let J be a periodic self-adjoint 5-diagonal matrix 
and 1Z be its spectral surface such that 1Z+ ~ C \ E. Since w{\) is single- 
valued in the domain, (2.7) and (2.9) imply that ujj = uj\ {Ej) + uJoo(Ej) 
are rational. Futher, due to (2.3) the function wz k is regular in the domain, 
moreover 1/ '(wz k )(Xo) and l/(wz k )(cc) are roots of the quadratic equation 

nfe-l ~ Tfk-l I ~ T-r2fc-l 

Thus (wz k )(X ) > and (wz k )(co) > 0. Since w = t(b Xo boo) k , t G T, the 
ratio (zb\ boo) k (oo) / (zbx b oo ) k (Ao) is also positive. That is e 2mk ^ = 1. And 
this finishes the necessity part of the theorem. 
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In the opposite direction, for the given system of cuts we define J 
according to Lemma 2.3. It remains to check that J is periodic. Let 
w = t^(b\ b oo ) k . Since n k (j) = 1, for all 7 € 7Ti(C \ E), it is single- valued. 
Note that w is normalized by the condition (wz k )((X)) = \(b\ b oo z) k (oo)\ > 0. 
We claim that 

we n = e n+2 k- (2.13) 

For odd n (2.13) holds automatically. For even n we should take into account 
that in addition t^t\ o = e~ 2mk ^ = 1. Thus (2.13) defines the shift operator. 
Since the multiplication operators by z and w commute, we have JS 2k = 
S 2k J. Therefore J is periodic. □ 

Now, let us restrict ourselves to the real case, i.e., C2 = c\ or both critical 
values are real. Without lost of generality C2 = c~\ = 2i or C2 = — c\ = 2. 
Thus, according to (2.5), 

z = \-\ (2.14) 
A 

in the first case, and 

z = \ + \ (2.15) 
A 

in the second one. 

In the case (2.14), z _1 (R) = R, thus E is a system of intervals on the real 
axis. Since z _1 (oo) = {0, 00}, E is subject to the restriction {0, 00} <f_ E. 

If z is of the form (2.15), then z _1 (IR) =1UT, that is, E is a union of 
real intervals and arcs of the unit circle, and again {0, 00} <£. E. This case 
under the additional assumption E CT leads to periodic CMV matrices [6] . 
Indeed, in the current case the multiplication by A is also well defined and 
represents a unitary matrix A such that 

J = A + A' 1 = A + A*. (2.16) 

As it was mentioned this functional model is the same as that related to 
periodic and almost periodic CMV matrices, see e.g. [4]. Conversely, having 
a periodic CMV matrix A we obtain the periodic self-adjoint J of the class 
by (2.16). 

Similarly, in the case (2.14) the multiplication by A leads to the self- 
adjoint operator A such that 

J = A-A- 1 , (2.17) 

where A -1 exists and corresponds to the multiplication by 1/A, i.e., to a pe- 
riodic SMP matrix. Further details of the corresponding functional model 
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are discussed in Section 4. Note that the case (2.15) under the additional 
assumption £cl leads to essentially the same class of self- adjoint opera- 
tors. 

3 Proof of the Main Theorem 

Let E = [bo, ao] \ U^ =1 (aj, bj) be a system of intervals on R, recall {0, oo} € 
C \ E, z = A — 1/A. We apply Theorem 2.1 in the current case. As it is 
well known the Green function (say with respect to infinity) is represented 
by the hyperelliptic integral, see e.g. [8, 7], 



where M K+ \ is a monic polynomial of degree x + 1. Note that the residue 
of the corresponding differential at the origin is —1. 

3.1 Spectrum of SMP matrices for the Stieltjes class 

Let us consider the simplest case E C M+ or E C M_ (we can say that 
the spectrum is on the upper (lower) sheet of 1Z + ). The strong Stieltjes 
moment problem is related to measures supported on the positive half-axis 
[3]. The shape of the sum G(X, oo) + G(X, 0) on M \ E is shown in Fig. 5. 
It implies immediately that all the zeros of the polynomial M x+ i in (3.1) 
are real in this case. Indeed, each gap (aj,bj), j > 1, contains at least one 
critical point; there is a critical point between — oo and 0; the total number 
of critical points is k + 1. Therefore 



is the Schwarz-Christoffel integral, which maps conformally the upper half- 
plane HI onto the (generalized) polygon in Fig. 6. 

According to (2.7) w(X) = te lke ( x \ t g T. Let Cbj be the coordinates of 
the base of the slits. Then w (A) is single valued in C \ E if and only if ujjk G 




Therefore for the sum of the Green functions we have 




(3.1) 




(3.2) 
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Figure 5: Go + Goo on R: the spectrum is on the upper sheet 
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Figure 6: Image of the Abelian integral 
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Figure 7: A complex critical point 



7rZ for all j. It remains to mention that due to the chosen normalization 
for b\ and the product bx (X)b oo (X)(X — 1/A) is positive at infinity and 
negative at the origin, that is £ = 1/2. Thus we can parametrize the spectral 
sets of periodic SMP matrices in this case by sufficiently simple domains 
shown in Fig. 6 with rational cu/s (quite similar to the Jacobi and CMV 
matrices cases). 

3.2 Complex critical points and three real critical points in 
the same gap 

The situation changes dramatically as soon as £ (aj,bj), j > 1, that is, 
E = E- U E + , E± C M±. Still all gaps, except for (aj,bj), should contain 
a critical point, thus M x+ \ has at least x — 1 real critical points. However, 
the positions of two remaining critical points are not a priory fixed. 

First, we consider the case when two remaining critical points are 
complex (io and JIq, Im/xo > 0. Let us consider 0(A) in the upper half- 
plane H. Since locally 6(X) = 9{hq) + C(A — ^q) 2 + . . . , there exist two 
orthogonal directions where Re d6 = 0. Moreover for one of them Im 6 has a 
local minimum at hq and a local maximum for another one. We define the 
curve 7, fio G 7, by the condition Re d6 = 0, such that Im 8 increases. Since 
there is no other critical point in EI this curve should terminate on the real 
axis. Note that Im 9(\) decreases as A approaches E. If so, in the gaps 7 may 
approach either a critical point or and 00. The first case is also not possible 
since the critical point is a local minimum for Im 6 = G\ + Goo along the 
real axis, thus it should be local maximum in the orthogonal direction 7. 
But along 7 it increases. Thus, 7 terminates at and 00, see Fig. 7. 

As the result we get HI \ 7 = H_ U H + such that E± C dW±. Let 6(\) = 
k9(X). Inspecting the boundary behavior of the given analytic function we 
obtain that it maps conformally W± onto D± shown in Fig. 8, where the 
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Figure 8: D± = 0(M±) regions for a complex critical point 

point —wo + itiQ corresponds to the critical point no. 
Now we define 



U (A) = y-i6(\) - icoo - ^ for A G H+ (3.3) 
Here we assume that Imu(A) > 0. Similarly we define 



u(\) = -yj-i0(\)-iuo-h% forAeH_ (3.4) 

and in this case Imu(A) < 0. In this way we get the regions . Since 

9 = -(2RenImu + oj ) +i((Reu) 2 - (Imn) 2 + hi), 

these regions are bounded by hyperbolic curves (1.1)-(1.4), see Fig. 1. Glu- 
ing the images along the curve 7 we obtain the conformal mapping of the 
upper half-plane H onto the special comb domain II = II + U II_ U R. 

Conversely, for the region II described by these equations we define a 
conformal map u : H — >■ II , it(0) = +00, n(oo) = —00, and set 

w (\) = e -(« 2 (A)+^+M.) 5 z = \-l (3.5) 

A 

Then, the set E corresponds to \w\ = 1. Since the base of the slits for D± 
are of the form ir£, w extended in the lower half-plane is single-valued in 
C \ E. Finally, w is real for A G R \ E, that is, ^ is rational. Based on 
Theorem 2.1 we conclude that this domain can be associated to a periodic 
SMP matrix. 

Let us turn to the case of three real critical points in the same 

gap. In this case H can be decomposed into three pieces. Let fi± < /xq < fJ-2 
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Figure 9: Three real critical points in the same gap, H + = Ho+ U 
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Figure 10: D± = 9(M±) three real critical points in the same gap 



be critical points in the gap (dj, hi). Note that with necessity fii and //2 are 
points of local maximum and Im# assumes a local minimum at /j,q in this 
interval. Therefore there are directions 71, 72 orthogonal to the real axis at 
fii and H2 respectively such that Im9 increases. Arguments like the above 
show that these curves, HedQ = 0, terminate at and 00, see Fig. 9 

In each of them, 9(A) represents a conformal mapping, see Fig. 10. In 
this picture — ujq + %h\, —ujq + ihf, and — uj$ + ih\ are images of the critical 
points no, fii, and ^2 respectively and wo = tt£o- We make the change of 
variable (3.3), (3.4), having in mind that now H+ or HL consists of two 
components. We arrive at the parametrization of the spectral curve by the 
domains of the form Fig. 2 such that 



h 



(4>) 



h 



h 2 



in (1.7). 
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As before, starting from a region II, by (3.5) we arrive at the set E and 
the domain C \ E c± 1Z + which corresponds to a periodic SMP matrix. 



3.3 Other cases 

In the previous subsection we considered critical values in two main generic 
positions. Now let us list the remaining special cases. 

1. For two complex critical values, if ujq = tt£q then one of the cuts in 
(1.3) and one in (1.4) degenerates to the intervals on the imaginary axis. 
The length of such a cut can not be arbitrary long, thus ti£ and hk + £ are 
subject for the conditions (1.5). As soon as one of these values approaches ho 
two complex critical values, from the upper and lower half-planes, approach 
the critical value in the corresponding gap. In the limit we have the critical 
value of multiplicity 3. The same special case can be obtained when two 
critical values [i\ and ^2 tend to /io, correspondingly h± — > 0, see Fig. 2. 

2. The case of a critical point of multiplicity two and a simple critical 
point in a gap corresponds to hf {,) = 0, h^ o) > or h (A)) = 0, h^ o) > 0. 

3. Two critical points (or one critical point of multiplicity two) may 
appear in the interval which contains zero or infinity. The domain II looks 
similar to that one shown in Fig. 2, but the degenerated hyperbola corre- 
sponds to the most left (or right) position, i.e., £q = or £q = m. 

4. It was assumed that m < k. If m > k the domain II in Fig. 1 remains 
the same, but we switch the normalization conditions to u(0) = —00 and 
u(oo) = +00. 

5. In the Stieltjes case, subsection 3.1, the spectral curve was described 
by a simpler domain, Fig. 6. By (3.3) it can be transformed to a II region 
bounded from below by the real axis. 



4 Functional model for periodic SMP matrices 

Let C \ E correspond to a periodic SMP matrix. We define 

For the reader's convenience we prove here a known lemma, see e.g. [4]. 
Lemma 4.1. The following identities hold true 

C(a) := Vl-W«)l 2 = MO) j^ (o y J (4-2) 
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and 



Proof. A{a) is defined by the following orthogonal decompositions 

fcoiC^ = ^/l-\A{a)\*K% Ajajb^K^ 

Indeed, .4(a) = and we get (4.2). 

Since, in addition, 

tf£ - AjajKg + ^pHPboC " 1 

600^0^ = \A-|-4(«)l 2 ^o 4(a)& iC° 1 
we have (4.3). 



□ 



In what follows without loss of generality we assume that = 1 in 
(2.12). In the given case Ao = 0, so to is the new notation for t\ , and this 
is not the same as the initial t( y Since 

(feooMM = (6oo A) (oo)6 (oo), (&ooM(0) = -(&o/A)(0)6oo(0), 

we have too = </>oo/|0oo| and to = -<po/\<t>o\, where 

/KA \oo), ^0= (A)(0)- 



bo J V^ooA, 
Also, recall that too I to = e 2m ^. 

Theorem 4.2. The multiplication operator by A with respect to the basis 
(2.10) is a periodic SMP matrix A = A(a, t(o)) with the following coefficients 

P2n =0oot(o ) e 2 ^ n ^(a^- ri )S~ 1 (a/^" n )C(a^oo / u- n )S(a A i 00/U - n ) 

P2n+i = - <Aoot(0)iooe 2 ^ ri C(a^- n )S(a Ai - ri )^(a^ ooA i- n )S- 1 (a^ oo ^- ri ) 

Q2n = ~ </>oo^(a / ti" n )S" 1 (a^ _Tt )^l(a^ 00 ^- n )S(a^oo/U~ n ) 

(4.4) 

and 

vr 2n+ i =t (0 )t 00 e 2 ^>oC(a Ai - n )B(a Ai -")^(a^ 00 ^-( n+1 ))^- 1 (a^ 00/ u-( n+1 )) 
vf 2 n+2 = - t(o)e 2 ^ (n+1 Vo^(a^ n )S- 1 (a^-")C(a / u 00 /i-(" +1 ))fi(^ 00/ u-^ +1 )) 
a 2 n+i = - ^o^(aM~ n )^~ 1 («Ai~")^(«//ooM- (n+1) )^(«Moo/x- ( " +1) ) 

(4.5) 
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Proof. Note that b^XU^ 6 i? 2 (a / u 00 ) and it is orthogonal to the subspace 
bob^H 2 (afiQ 1 fi^) . Therefore, in fact, we have the three-terms recurrence 
relation, 

Ae2n = p2n^2n-l + 02n^2n + P2n+ie2n+l • 



Moreover 



P2n = t(o ) e 2 ^"(A6 00 )(oo) K ^ {CO) 



and 



= {XK7 - _ (A&oo)(oo)irg^ " C -(oo),C") 

&oo#S^ (OO) 

(Ab 00 )(oo) irff°°"~"(0) K^ n (oc) 
P2n + i =<At (0) to C'~ 1 M&' i ° V "> 



t (0) e lW f„o ^^-(oc) } 



^ (A600K00) C-^jO) 6 (oo)^r° M (00) 
^ (0) °° ^(0) X«~"(oo) iC"» 



In its turn, 

1 



where 



^2n+l — K2n+\Z2n + (^2n+1^2n+l + 7T2 n + 2 e 2n+ 2 , 



Vf2n+1 = t( )tooe- 2 ^" (0) 



a 2n+1 =(^sr° v " n ,ifsr°" v " n ) 



" ( A^ "UJ (0) 6o^-"(0) 00 ( ) ' 00 } 

Voc) ^-"(0) ^oV" (oo) ' 



-l..-r 
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and 

77 ■') . . : — ( ^ 



„ /^-n-lr/% 1 f " + +-n-li j>-«M V ™\ 

7T2n+2 = \T*oo A oo ; *(0)*0 °oo^ / 



- _f -27ri£(n+l) /'WooVm 1 (00) (0) 

(0)e v a ; w moo) ^- ( o )lC o-v- (00) - 

Now by making use of (4.1), (4.3), (4.2), we obtain (4.4) 
<?2n = - </>oo^(a/^oo^")^(aAtoo^~")^(a^~ n )^~ 1 (a/ u ~ n ) ) 

^ _J, + ^wi^n, fn \ (0°) 

p 2 n-0oot (O) e M0)^_„ M 

() (oo)iC^ (00) 1 

=0oot(o)e 27r ^(a Ai -")fi- 1 (a/i- n )C(a^ oo/[ /-")S(a /Uoo /i- n ), 



tfST" (00) ff?"~"(0) 

(-(0)*oo 

as well as (4.5) 



t(0)tooe- 27r ^(a^ oo/U -™)S(Q^ oo ^ n )C(a^ ri )fi- 1 (aAi- n ) 



7T 2 n+i = i(o)*ooe ? </>o ^—^ 00(00) — — — — — 

KZ° M (00) (00) (0) 

= t ( o)iooe- 2 ^ n 0o^(o^H^6(^oo^ 



cT2n+i = -0o^(a^ ri )S- 1 (a^ n )^(a / u ooA i-("+ 1 ))S(a^ ooM -( n+1 )) 
W2 = -t (0) e 2 ^( n+1 Vo^(a/i~ n )^- 1 (a/"" n ) 
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-t (0) e 2 ^( n+1 )^(a /[ i- n )fi- 1 (a Ai - n )C(a/i oo ^(" +1 ))S(a Aloo ^( n+1 )) 



□ 



Remark 4.3. The structure of the reproducing kernels on the hyperelliptic 
Riemann surfaces is well known, see e.g. [7]. In particular, indeed -fT^O) = 
0, i.e., A(a) = 0, for some a. According to (4.4) and (4.5) it means that the 
corresponding A = A(a, i( )) may degenerate, that is, q2 n or e>"2n-i vanishes 
for some n. Nevertheless all entries of A and A^ 1 have perfect sense. For 
example, 

P2nP2n+l 
T2n+1 = 

=|0 oo |C(a/i~ n )S(a Ai - n )C(a^ oo/ u- n )S- 1 (a^ ooM - n ), 

^2n-1^2n 

~P2n = 

— C2n-1 

=|0o|C(a^ ( ^ 1) )^(^~ {ri ~ 1) )C(a/ioo^")6" 1 (^oo^- n ), 

where 

Ae2n-\ =^2n-ie2n-3 + p2n-l&2n-2 + <?2n-ie2rt-l + P2v,e2n + ^2n+ie2n+l, 
A~ 1 e2n =P2n^2n-2 + ^2n^2n-l + ^2n^2n + ^2n+1^2n+l + p2n+2&2n+2- 
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